Glasslike vs. crystalline thermal conductivity in carrier-tuned Ba 8 Gai 6 X 30 clathrates 

(X = Ge, Sn) 
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The present controversy over the origin of glasslike thermal conductivity observed in certain 
crystalline materials is addressed by studies on single-crystal x-ray diffraction, thermal conductivity 
k(T) and specific heat C P (T) of carrier-tuned BasGai6X3o (X = Ge, Sn) clathrates. These crystals 
show radically different low-temperature k(T) behaviors depending on whether their charge carriers 
are electrons or holes, displaying the usual crystalline peak in the former case and an anomalous 
glasslike plateau in the latter. In contrast, C P (T) above 4 K and the general structural properties are 
essentially insensitive to carrier tuning. We analyze these combined results within the framework of 
a Tunneling/Resonant /Rayleigh scatterings model, and conclude that the evolution from crystalline 
to glasslike k(T) is accompanied by an increase both in the effective density of tunnelling states and 
in the resonant scattering level, while neither one of these contributions can solely account for the 
observed changes in the full temperature range. This suggests that the most relevant factor which 
determines crystalline or glasslike behavior is the coupling strength between the guest vibrational 
modes and the frameworks with different charge carriers. 
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I. INTRODUCTION 

The understanding of thermal conductivity behavior 
k(T) is of direct interest to any research involving the 
discovery, design and development of materials for ther- 
moelectric conversion applications, where k(T) should 
be as small as possible, while at the same time ther- 
mopower S(T) and electrical conductivity cr(T) should 
be large. In the semiclassical theory for electron and 
phonon transport several mechanisms are known as 
contributors to heat conduction/phonon scattering in a 
material, consequently affecting its overall thermal con- 
ductivity. 

In metals, heat conduction by charge carriers is 
the largest contribution, and is well described by the 
Wiedemann-Franz law k c = LaT, which directly relates 
the carrier thermal conductivity k c with an appropriate 
Lorentz number L~2-3x 1CT 8 WQ/K 2 , the electrical 
conductivity a and the temperature T. Due to their typ- 
ically large charge carrier densities n c , metals have large 
a(n c ,T) and thus large k c in the range of 50 — 500 W/m K 
at room temperature. 

Conversely, semi-metallic, semi-conducting and insu- 
lating compounds have low n c and <t(T), therefore small 
and often negligible k c (T) and the overall heat conduc- 
tion behavior k(T) near room temperature is in the 
range of 10 — 50 W/m K, governed mostly by con- 
tributions kl(T) arising from the crystal lattice. At 
T — > 0, k(T) — > from basic thermodynamic principles, 
so within the first few Kelvins k(T) increases quickly as 
a power law T r , with 1 < r < 3 depending on which 
phonon scattering mechanisms dominate at low temper- 
atures. At higher temperatures the phonon scattering 
is generally described as governed mostly by umklapp 



processes, for which the Debye approximation approach 
shows a decrease with a T -1 dependence. Therefore, at 
some intermediate temperature usually around 10 — 50 K, 
a characteristic "crystalline peak" is observed in k(T) 
for common compounds, due to the crossover from one 
regime to another. 2 

The peak usually appears equally in polycrystalline 
materials since grain boundary scattering is in general 
a minor contribution^^ unless the average grain size be- 
comes very small or the temperature very low. How- 
ever, glasses are an exception to the above because of 
two basic factors: the very low mean free path for both 
electrons and phonons, and the presence of low-energy 
tunnelling states (TS), i.e., different localized potential 
minima for atomic positions in their amorphous distribu- 
tion of nuclei^ This class of materials shows extremely 
low heat conduction and a universal behavior of k(T): 
the lowest temperature behavior (up to ~ 1 K) rises as 
T 2 due to scattering by the tunnelling states, then levels 
off as a characteristic intermediate temperature plateau 
(attributed to Rayleigh scattering). Above ~ 20 K it re- 
sumes a slow increase, until roughly levelling off again at 
higher temperatures. 

In more recent years, the search for new and potentially 
useful thermoelectric materials^ has led to the discovery 
of compounds that not only have unusually low thermal 
conductivity, but whose general behavior resembles that 
of a glassy material despite the fact that they are true 
(albeit disordered) crystalline lattices^ A prominent ex- 
ample is the intermetallic compound SrgGai6Ge3o, with a 
filled type-I clathrate structure^ (6 larger X24 cages form- 
ing tetrakaidecahedra plus 2 smaller X20 cages forming 
dodecahedra, see Fig. ^) for which k(T) was first mea- 
sured by Nolas et al. in 1998. 3 A model was proposed^ 
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FIG. 1: (Color online) The two cages of the type-I clathrate 
structure adopted by BasGaieGeso- If we consider the cage 
"construction unit" as 4 atoms connected in zig-zag from top 
to bottom, then the larger X24 cage (left) is made of 6 such 
units and the smaller X20 cage (right) is made of 5 units. 



to explain this material's glasslike behavior, based on the 
idea that TS exist for the Sr(2) guest ion in the X 2 4 cage^ 
to which it is rather loosely bound because of an ion- 
to-cage size mismatch. A combination of phonon scat- 
tering by TS, resonant scattering on large, Einstein-like 
localized vibration modes (guest rattling) and Rayleigh 
scattering was used to adequately reproduce the exper- 
imental k(T) behavior (henceforth we will refer to this 
combination as the TRR model). Later investigations 
amply demonstrated a splitting of the Sr(2) site into four 
off-center positions jiiLlLH among which the ions could 
indeed tunnel. 

As other clathrate compounds started being investi- 
gated, the TRR model was challenged by at least two 
other models. One proposes that the tunnelling states 
are not required, only an off-center vibration of the guest 
ionsjl&U and another proposes that the guest ions don't 
play a major role at all at low temperatures, but rather 
it is the phonon scattering on charge carriers that leads 
to the glasslike behavior ^iLiLiS 

In this work we address the issue by performing sin- 
gle crystal x-ray diffraction (SCXRD), thermal conduc- 
tivity k(T) and heat capacity C P (T) experiments on 
BagGaieSnso (BGS) and BasGai6Ge 3 o (BGG) crystals, 
which have been tuned through the crystal growth pro- 
cess to display n-type or p-type majority charge carriers 
as a result of small imbalances in their Ga:Ge or Ga:Sn 
ratios>i2i2£ By analyzing the differences and similarities 
between the behaviors of these samples, we can test the 
applicability of the various models proposed to explain 
the origin of unusual glasslike behavior, in this case ob- 
served for p-type samples whereas the n-type samples 
show the normal crystalline peak. 




FIG. 2: (Color online) Irregular cage of the type-VIII 
clathrate structure adopted by BasGaieSnso. The four 
smaller cage spheres represent the 8c site, preferentially oc- 
cupied by Ga atoms. 



TABLE I: Average Ba:Ga:X content (X = Ge, Sn) in the 
four measured crystals as determined by electron-probe mi- 
croanalysis. 



Sample 


Starting Flux 


Crystal 


Name 


Composition 


Composition 


n-BGS 


8 : 16 : 60 


8.0 : 15.98 : 30.02 


n-BGG 


8 : 38 : 34 


8.0 : 15.94 : 30.06 


p-BGS 


8 : 38 : 30 


8.0 : 16.14 : 29.86 


p-BGG 


8 : 38 : 30 


8.0 : 16.10 : 29.90 



II. EXPERIMENTAL DETAILS 

Growth details of the batches used in this work are de- 
scribed in previous papers Single crystalline poly- 
hedrons of 3-10 mm in diameter were obtained by a self- 
flux method. The carrier type is tuned by choosing Ga or 
Sn flusiS in the case of BGS or by adjusting the relative 
Ge content in the initial mix with Ga fluj*2£ in the case 
of BGG. The batch name, flux composition and crystal 
composition determined by a JEOL JXA-8200 electron- 
probe microanalyzer (EPMA) are summarized in table [I] 
The composition values are averages over 10 regions of 
each crystal, although there are random fluctuations of 
up to ±0.1 throughout the crystals. The values shown 
for the p-BGG sample should be considered a correction 
to those published in refs. |2(] and |2^, since this was a 
more careful evaluation made on the same batch. As ex- 
pected from charge balance principles, Ga-rich samples 
show p-type carriers while Ge-rich samples show n-type 
carriers. 

Thermal conductivity experiments were performed us- 
ing a steady-state method on home-made systems, in 
the range of 0.3-300 K (BGG) and 4-300 K (BGS), al- 
though reliable data is only obtainable up to about 150 K. 
At higher temperatures, thermal losses by radiation and 
wire conduction prevent the correct measurement of the 
intrinsic sample properties. The electronic contribution 
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FIG. 3: (Color online) Temperature dependence of the ther- 
mal conductivity of BasGaieGeso and BasGaigSnao with dif- 
ferent carrier types. Solid lines are best fits of the TRR model 
as described in the Discussion section. 



k c (T) of all samples estimated by the Wiedeman- Franz 
law is negligible up to 100 K, so the measured k(T) is 
equated to the lattice contribution kl{T). Heat capac- 
ity was measured using a Quantum Design PPMS with 
its standard thermal-relaxation method, in the range 
0.4 < T < 300 K. 

For SCXRD experiments, broken pieces of n-BGS 
and p-BGS with approximate dimensions of 0.1 x 
0.1 x 0.1 mm 3 were selected. The diffraction data 
were collected with a Rigaku R-AXIS RAPID imaging 
plate area detector using graphite-monochromated Mo 
K a radiation. Refinements were performed using the 
CrystalStructure 23 software. Structures were solved by 
direct methods and expanded using Fourier techniques. 
All sites were assumed to be fully occupied. 



III. EXPERIMENTAL RESULTS 

A. Thermal Conductivity 

The lattice thermal conductivity kl (T) of all four sam- 
ples is shown in Fig. [31 (symbols are the as-measured ex- 
perimental data points and solid lines represent fits to 
the data by the TRR model which will be detailed in the 
Discussion section). At 100 K and above (not shown), 
k l (T) for BGS is roughly half that of BGG (about 1 and 
2 W/m K respectively). This can be understood as a con- 
sequence of three main factors: (i) if the rattling of the 
guest ions is the main contributor to the unusually high 
phonon scattering level in these materials^ the larger 
cage size in BGS leads to larger rattling of the guest 
ions; (ii) in the BGS unit cell all 8 guest ions vibrate 
with equal intensity (single crystallographic site for Ba 
in the Type- VIII clathrate structure) , while in BGG only 
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FIG. 4: (Color online) Heat capacity C p of BagGaieSn3o 
and BagGai6Sn3o with different carrier types, presented as 
Cp/T 3 vs.T. Solid lines are best fits of the Einstein model as 
described in the Discussion section. 



the 6 guest ions inside the X24 cages show large rattling; 
and (iii) the heavier Sn atoms produce lower frequency 
phonons, which are more easily scattered. 

Below 100 K, each sample behaves quite differently, de- 
pending more on the carrier type than on the compound. 
The two n-type samples increase towards a peak, while 
the p-type samples remain at a plateau, smaller by a fac- 
tor of 3 — 4 in value than the n-type counterparts near 
the peak. Below 10 K, kl(T) for all samples decreases 
fast and, in the case of BGG which was measured to 
lower temperatures, a gradual crossover to a T 2 regime 
is clearly observed. This implies a phonon mean free 
path inversely proportional to frequency^ which is the ex- 
pected dependence when phonon scattering by tunnelling 
states is dominant. The T 2 behavior contrasts with a 
previously reported result showing a T dependence for 
p-type BGGiiSi 



B. Heat Capacity 

The data points in Figs. 0k and show the as- 
measured specific heats C P (T) for the BGS and BGG 
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TABLE II: Summary of crystallographic parameters from the structural refinement of a n-type BagGaigSnao single crystal. 
Space group I43n (No. 217), a = 11.586(1) A, Z = 1, R = 0.009, R w = 0.009, B eq = (8tt 2 /3) J2i J2j U H < a>iOj. 



Atom 


site 


X 


y 


z 


B eq (A 2 ) 


occupancy 


Ba(l) 


8c 


0.68490(5) 


0.31510(5) 


0.31510(5) 


3.32(1) 


1 


Ga(l) / Sn(l) 


12d 


0.5000 


0.0000 


0.2500 


1.65(1) 


0.184(12) / 0.816(12) 


Ga(2) / Sn(2) 


2a 


0.5000 


0.5000 


0.5000 


0.97(2) 


0.158(22) / 0.842(22) 


Ga(3) / Sn(3) 


24g 


0.41549(4) 


0.14887(4) 


0.41549(4) 


1.333(9) 


0.314(8) / 0.686(8) 


Ga(4) / Sn(4) 


8c 


0.36558(6) 


0.36558(6) 


0.36558(6) 


1.07(2) 


0.766(12) / 0.234(12) 



TABLE III: Summary of crystallographic parameters from the structural refinement of a p-type BagGaigSnso single crystal. 
Space group I43n (No. 217), a = 11.587(1) A, Z = 1, R= 0.0218, R w = 0.0157, B eq = (8tt 2 /3) £\ £\ U lj a*a*a l a j . 



Atom 


site 


X 


y 


z 


B eq (A 2 ) 


occupancy 


Ba(l) 


8c 


0.68507(8) 


0.31493(8) 


0.31493(8) 


3.25(2) 


1 


Ga(l) / Sn(l) 


12d 


0.5000 


0.0000 


0.2500 


1.69(2) 


0.152(16) / 0.848(16) 


Ga(2) / Sn(2) 


2a 


0.5000 


0.5000 


0.5000 


0.90(3) 


0.233(26) / 0.767(26) 


Ga(3) / Sn(3) 


24g 


0.41565(4) 


0.14836(4) 


0.41565(4) 


1.308(14) 


0.318(10) / 0.682(10) 


Ga(4) / Sn(4) 


8c 


0.36577(6) 


0.36577(6) 


0.36577(6) 


1.077(14) 


0.707(12) / 0.293(12) 



samples respectively, plotted as Cp/T 3 vs. T. This plot- 
ting style emphasizes the contributions of localized vibra- 
tions of guest atoms (Einstein oscillators), which appear 
as pronounced peaks over a "background" contribution of 
a Debye solid. For these samples, the charge carrier con- 
tribution is negligible above ~ 4 K, but responsible for 
the T~ 2 upward curvature upon cooling below this tem- 
perature. A more traditional plot of Cp/T vs. T 2 below 
4 K (not shown) is used to estimate with good accuracy 
the Sommerfeld coefficient 7 of the charge carriers and 
the Debye temperature Qd of the 46 framework atoms, 
and then subtract the host contribution Ct (dotted lines) 
in order to isolate the Einstein-like contribution C g of the 
guest ions (dash-dotted peak). 

Contrary to the heat transport data, in both BGS 
and BGG the heat capacity data show the same behav- 
ior above ~ 4 K for different carrier types. This result 
demonstrates there is no fundamental change in the en- 
tropic properties of these compounds within the range 
of deviations from stoichiometry studied. If the rattling 
behavior of the guest ion is not significantly changed for 
different carrier types in the framework, then it should 
be the coupling between the guest vibration and the 
frameworks with different carriers that changes. In other 
words, frameworks with holes have their phonon modes 
more effectively scattered by the Ba vibration than those 
with electrons. 



tling) of the guest ion in the oversized cage, however, no 
relevant difference is observed between the crystals. 

The resulting sets of data do not allow a detailed com- 
position analysis for comparison with EPMA results, be- 
cause the R factor was insensitive to occupation probabil- 
ity within deviations of ±0.2 from stoichiometry, but the 
framework sites show consistent preferential occupations 
for Sn(l) and Sn(2) in the respective 12c? and 2a crystallo- 
graphic sites, while Ga(4) has the preferential occupation 
of the 8c site (in accordance with the idea that the smaller 
Ga atom should more easily occupy the site with smaller 
bond distances between neighbors) and the 24g site re- 
mains more randomly occupied by Sn(3) and Ga(3). This 
is true for crystals with both types of carriers, the only 
consistent and relevant difference we could find in these 
refinements was a larger relative occupation of the 2a site 
by Ga(2) for the p-BGS samples (the top atom in Fig.|2J. 
This could be where the "extra" Ga ions prefer to enter 
in Ga-rich samples, but whether or not this can have any 
influence on the overall guest/framework coupling would 
require more detailed investigation. 

IV. DISCUSSION 

We now present and discuss the models used to analyze 
the data in Figs. |21and2| The specific heat is expressed 
as a sum of 3 main contributions: 



C. Single-Crystal X-Ray Diffraction C p = C C + C D + C E , (I) 

. . . . where C c = 7T is the electronic specific heat of the 

Tables |II| and |III| summarize the refinement results charge carriers 
made for room temperature SCXRD data of n-BGS and 
p-BGS respectively. The anomalously large isotropic 

thermal parameter B eq of the Ba site compared to the „ 12ir 4 NDkB f &D / T 



Ga/Sn sites is a signature of the enhanced vibration (rat- D 5 J Q (e x — l) 2 
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with x = hck/ksT is the Debye model for the lattice 
specific heat of N D Debye oscillators per unit cell, whose 
numerical solution can be found in Solid State Physics 
textbooks, and 

is the Einstein specific heat of the i-th vibrational 
mode of any existing rattling ions. For our analysis we 
assume that the 8 rattling guest ions are sufficiently de- 
coupled from the 46 rigid framework atoms, so that we 
can make the association C g = Cei and C'h — C c + Co 
respectively. 

The solid line in Fig.^J, indicates the best fit of eq.|3]to 
the isolated Einstein contribution in p-BGS. It is impor- 
tant to emphasize that in this analysis for the type- VIII 
structure, the dimensionality and the number of Einstein 
oscillators are fixed at p — 3 and Ne = 8 so there is a 
single fitting parameter Qe for BGS, which alone governs 
all the peak characteristics - position, height and width. 
The fact that the best fitted curve with 0^ = 49.9 K so 
closely reproduces the experimental peak in all these as- 
pects is a solid testimony to how successful the Einstein 
model is in describing the vibrational behavior of the 8 
Ba ions in this compound. The data is actually slightly 
broadened with respect to the model, which may be the 
result of a narrow distribution of Qe around the mean 
value (a consequence of Ga/Sn site disorder), and/or a 
slightly anisotropic vibration of each guest ion in its re- 
spective irregular cage (resembling an ovoid with diam- 
eter varying between 7.3 — 8.2 A, see Fig. We will 
see next how anisotropic vibration plays a much more 
important role in the BGG compound. 

Contrary to BGS, the specific heat of BGG cannot be 
adequately fit with a single Qe and Ne = 8. If Ne is 
freed as a fitting parameter, the best fit naturally de- 
creases this to Ne — 6.1 (and Qe = 55 K), consistent 
with the fact that the 6 Ba(2) ions in the larger A24 
cages are the main rattlers, but the fitting quality is still 
not satisfactory. Until now, the usual approachi^2S to 
analyze the heat capacity of BGG has been to assign 
two different Einstein contributions (i = 2 in Eq. 
This results in excellent fits with Qei = 70 — 80 K and 
Qe2 = 30 — 40 K. However, the number of Einstein oscil- 
lators Nei and Ne2 results opposite to what one would 
expect if these numbers were to represent the two Ba 
sites, i.e., there is a greater number of Ba oscillators with 
larger O^i (Nei =6 — 9) than those with smaller Qe2 
(Ne2 = 1-5 — 2.0). This is difficult to physically justify, 
since larger rattling implies smaller Qe- 

We offer an analysis which better reconciles with the 
guest ions' known physical properties. The starting point 
is that, due to the X 2 4 cage shape (Fig. ^ left), the 
Ba(2) ions show a strongly anisotropic vibration with 
greater amplitude within the plane parallel to the cage's 
two hexagons^ Because the dimensionality p plays a role 



in the Einstein model (see Eq. |3J, at least two vibra- 
tional modes should be required to describe the Ba(2) 

site alone: in-plane {Q^ E2 j and out-of-plane (Q E2 ) re ~ 
spectively. In addition, a third vibrational mode (0_bi) 
is required to account for the smaller, but still Einstein- 
like, rattling of the two Ba(I) site ions in the X20 dodec- 
ahedra (Fig.^ right), which can be assumed isotropic. 11 
In this model, the dimensionalities and numbers of oscil- 
lators are predefined: p\Nei =3x2, p\n e2 = 2x6, 
P2^E2 = 1 x 6, so the fitting parameters are only the 
three Einstein temperatures, with the additional con- 
straint that Q E2 < Qe2i®ei- The best fitting of this 
model results in Q E i = 87.2 K, e|, 2 = 49.4 K and 
Q E2 = 87.1 K. The similarity between Q% 2 of BGG 
and G^ of BGS is reasonable, since the largest diame- 
ters of both cages are essentially the same (~ 8.2 A). 
The similarity between Q E2 and Qei in BGG is also 
reasonable since the X24 cage size in the out-of-plane 
direction is very close to the X20 cage size (~ 5.5 A). 
This means a further simplification can be made in the 
model by assuming only two parameters Qei and 0^2 
with p\Nei = P2NE2 = 12, where Qei represents the 3D 
vibration of the Ba(l) ions and the ID out-of-plane vibra- 
tion of the Ba(2) ions; while Qe2 represents the larger, 
2D in-plane vibration of the Ba(2) ions. This results in 
the solid curve shown in Fig. and, as with BGS, the 
data for BGG is only slightly broadened with respect to 
the model. 

With the heat capacity parameters determined, we 
now analyze the lattice thermal conductivity kl of all 
samples, using the same procedure applied previously for 
analysis of the n-BGS sample^ which is in turn based 
on the TRR model initially used in ref. @ to describe 
SrgGaigGe3o. In the semi-classical theory, kl is given by 

*l = -£J <L,[C L (u,T)vt\, (4) 

where Cl{u),T) is the phonon specific heat, lu e is the 
Debye frequency, v is the average sound velocity and I is 
the phonon mean free path, which must be averaged over 
all major contributing scattering mechanisms. Thus, in 
the TRR model it is written as 



' — i^TS + ^Res + ^Ray) + ^nin- (5) 

The low-energy excitations of the guest ions tunnelling 
between localized states scatter phonons as 

TS \k B J \2k B T J 2 \tiu BT 3 J ' 

- . (6) 

where A and B are microscopic parameters describing 
the tunnelling states characteristics^ At higher energies, 
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TABLE IV: Parameters used to generate the solid line curves 
in Figs. |3]and3] which best fit the respective experimental 
data set for lattice thermal conductivity and specific heat. 
See text for detailed descriptions. 



Symbol 


Unit 


n-BGG 


p-BGG 


n-BGS 


p-BGS 


A 


107(m K) 
1/K 2 


1.4 


2.5 


2.5 


17 


B 


0.1 


0.1 


0.1 


0.1 


A/B 


10 5 K/m 
l/(m s 2 K ) 


1.4 


2.5 


2.5 


17 


Ci 


0.2 


2.0 


0.7 


5.0 


Qei 


K 


87 


87 


50 


50 


Ti 




0.5 


1.5 


0.4 


1.5 


C'2 


l/(m s 2 K 2 ) 


0.2 


2.0 






&E2 


K 


49" 


49" 






r 2 




0.5 


1.5 






D 


K 4 /m 


0.85 


0.5 


2.8 


1.7 


7 


mJ/(mol K 2 ) 


6 


9 


1.3 


11 


e D 


K 


288 


288 


200 


200 


V 


m/s 


2898 


2898 


2250 


2250 



"Two-dimensional vibration (see text). 

phonons are scattered through a resonance effect against 
guest ion rattling as: 



r 1 

l Res 



E 



C,uj 2 T 2 



{ujf - uj 2 ) 1 + T^fto 2 



(7) 



where Ci and I\ are phenomenological parameters re- 
lated to a simple mechanical oscillator. 26 We also need to 
include the empirical but always present, frequency-only 
dependent Rayleigh scattering 



r 1 

L Ray 



D 



hw 4 



(8) 



and finally the last term l m i n — 1 A is the cut-off limit. 

Results from the best fits of the data shown in Figs. 
13 and 0| are summarized in Table IIVI The most rele- 
vant results in terms of comparing the p-type with n- 
type samples are the increase in the resonant scattering 
level (Ci), and in the TS scattering level. The latter can 
be expressed by the ratio A/B — h(hv) 2 /irks, which in 
glasses is essentially a measure of the subset density of 
tunnelling states n that are able to strongly couple to 
the phonons and effectively scatter them* 2 ^ Therefore, 
the increase in A/B observed upon changing from n-type 
to p-type cages does not necessarily mean the total den- 
sity of TS has increased, only that the existing states are 
more effectively coupled. 

An interesting exercise can be made to help understand 
the influence of these different contributions in the TRR 
model. If we begin with the fitting results for «l(T) of 
the n-type samples, it is impossible to fit the respective 
p-type Ki,(T) by increasing the intensity of only one of 
these contributions (TS or resonant). The TS are mainly 
responsible for decreasing the low-temperature Kl (T) up 
to the first few Kelvins, and by itself the TS contribution 



is incapable of changing the peak into a plateau. Con- 
versely, an increase in the resonant scattering level (based 
on fixed values of Ob = 50 K from heat capacity and in- 
creased phenomenological coupling strength parameters 
Ci), readily brings the peak down to a plateau/dip, but 
quickly loses its ability to follow the Kl{T) drop below 
about 10 K. Therefore, we may conclude that the TRR 
model adequately reproduces the entire range of kl(T) 
for all samples up to 100 K, provided that the coupling 
of the framework phonons with the guest ion tunnelling 
and rattling is increased in p-type samples. 

Let us now focus on some other proposals regard- 
ing candidate mechanisms for glasslike behavior in 
clathrates, which challenge the TRR model. First: is 
the presence of tunnelling states really necessary, or is 
off-center vibration a sufficient mechanism? The ques- 
tion was raised by Bridges and Downward^ primar- 
ily based on the existing data at the time, where Sr 
and Eu guests clearly show off-center sitesi2*ii*i£ and 
glasslike k(T), while Ba guests appear to show on-center 
vibration 10,11,12 (within experimental resolution) and a 
crystalline peak. Later studies demonstrated glasslike 
behavior for p-type BG G 15 i 20 and now for p-type BGS 
(this work), so this argument by itself is no longer valid, 
unless a closer look at the Ba vibration in these com- 
pounds through microscopic techniques shows that off- 
center vibration modes do exist for p-type samples (even 
if much smaller than for Sr and Eu guests) but not for 
n-type samples. Raman scattering and EXAFS studies 
are presently being conducted on our carrier-tuned BGG 
crystals, which may help clarify this issue. Still, good 
arguments were made by the authors in terms of describ- 
ing how off-center vibration does indeed help enhance the 
coupling between guest vibration modes and the frame- 
work phonons ^2*i4 

A second challenge to the TRR model is: can the shift 
from crystalline to glasslike behavior be explained solely 
by phonon scattering mechanisms within the framework, 
i.e., by interactions between phonons and charge carriers? 
This question was raised in a series of papers by Bentien 
et tilJ-LlLl— which we now discuss. 

The first worki^ called attention to an observed k(T) oc 
T 1,5 dependence at low-temperature for p-BGG and a 
kink in their data at about 2 K (neither of which were 
reproduced with our crystals). They also pointed out 
that the phonon-charge carrier mechanism could not ex- 
plain the lowering of kl(T) above ~ 15 K, so the reso- 
nant scattering on the guest vibration was once again in- 
voked, but to account for only the differences above 15 K. 
The second work 16 compared several polycrystalline sam- 
ples of type- VIII and type-I EusGai6Ge3o (a-EGG and 
/3-EGG respectively^ all with n-type carriers), clearly 
demonstrating that /3-EGG shows glasslike kl(T) while 
a-EGG does not. The difference was interpreted in terms 
of changes in the band structure, with a much enhanced 
effective mass m* found in /3-EGG. However, the cage 
sizes and shapes are also quite different between these 
two structures. The type-I A24 cages are essentially the 
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same size for all Ge clathrates (5.5 x 8.2 A, see Fig.^) but 
the type- VIII cage in a-EGG (6.7 x 7.5 A ovoid similar to 
Fig. [21 is significantly smaller than that of BGS, so any 
change in kl(T) can also be argued or modelled in terms 
of changes in the Eu vibration modes and their coupling 
to the framework. Unfortunately a-EGG samples with 
p-type carriers are as yet unavailable, but it wouldn't be 
surprising if they showed glasslike kl{T) as we found in 
p-BGS. The third and more recent worki& shows results 
for Ba8Ni r Ge46_ :r similar to what we have obtained here 
for BGG and BGS, therefore the same analysis and dis- 
cussion we have conducted here can also be applied to 
those results. 

Still, it is obvious that the influence of charge-carriers 
cannot be neglected with respect to their density n c , ef- 
fective mass m* , electronic mean free path l c , etc. It 
is quite clear from our measurements and all previously 
reported data on Ba-filled clathrates, that the p-type car- 
riers are playing a relevant role in producing an increased 
phonon scattering in these compounds, which we view as 
yet another additional factor capable of contributing to 
lower kl (T) , possibly through direct interaction with the 
phonons, but especially by mediating an enhanced cou- 
pling of these with the guest vibration modes. A few 
brief examples for such mediation possibilities are: 1) n- 
type frameworks could allow a greater degree of coherence 
in the vibrations of neighboring Ba guests than p-type 
frameworks, which would lead to larger mean free paths 
and less effective scattering; 2) Since the type and density 
of charge carriers result from stoichiometry imbalances, 
they may affect the framework rigidity at certain sites, 
and therefore how easily it can couple with the rattler 
ions. 



V. CONCLUSION 

We have succeeded in growing large single crystals of 
BasGai6Sn 30 and Ba 8 Gai6Ge3o with both n-type and p- 
type majority carriers, and found that these compounds 
show low temperature lattice thermal conductivity be- 
havior strongly dependent on the carrier type. A shift 
from crystalline to glasslike behavior is observed for both 
compounds when changing the majority carriers from n- 
type to p-type through composition tuning. These differ- 
ences can be mostly reproduced by an increase in reso- 
nant scattering, however, an increase in both resonant 
and tunnelling scattering levels are required to repro- 
duce the full set of data below 100 K. Heat capacity 



and single-crystal x-ray diffraction data indicated that 
these increases are not the result of any major change 
in the guest ions' vibrational behaviors, therefore a more 
effective coupling of the frameworks with p-type carriers 
to the TS and rattling vibrations of the guest ions is the 
most likely mechanism. The T 2 dependence in kt, (T) ob- 
tained at lowest temperatures for both n-type and p-type 
BasGai6Ge3o indicates that tunnelling states should be 
present for the Ba(2) ions in this compound, therefore its 
mere presence is insufficient to guarantee glasslike kl(T). 

In fact, our results indicate that the various proposed 
mechanisms which may lead to glasslike behavior are 
all partially correct and at the same time incomplete. 
The general scenario that we see emerging can indeed 
be expressed as: it's all about the coupling. For reasons 
that still need to be explained microscopically, the n- 
type frameworks are more weakly coupled to the guest 
vibration modes than the p-type frameworks. Thus, the 
Ba ions' smaller and (almost?) on-center vibration is 
not coupled strongly enough to the n-type framework 
phonons to produce the glasslike behavior, but the p- 
type framework crosses the necessary coupling strength 
threshold to achieve this scattering regime. In contrast, 
Sr and Eu ions in the type-I Ge clathrates have clearly off- 
center and larger rattling, capable of a strong enough cou- 
pling even with the n-type frameworks to produce glass- 
like behavior (no p-type frameworks have been reported 
yet for these compounds). In a series of carefully tuned 
Ba-based clathrates it should be possible to observe a 
continuous transition from glasslike to crystalline kl(T). 
Likewise, in a series of n-type (Sr,Eu)-based clathrates 
the same continuous transition should be observed not 
from carrier tuning, but from a physical or chemical re- 
duction of cage size to dampen the off-center vibration 
level. 
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